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Abstract 

A spin filtering device using quantum spin interference is theoretically proposed in a 
GaAs/AlGaAs electron gas that has both Rashba and Dresselhaus spin-orbit couplings. The device 
achieves polarized electron currents by separating spin up and spin down components without a 
magnetic field gradient. We find two broad spin filtering regimes, one where the interferometer has 
symmetrical arms, where a small magnetic fiux is needed to achieve spin separation, and the other 
with asymmetric arms where the change in path length renders an extra phase emulating the effects 
of a magnetic field. We identify operating points for the device where optimal electron polarization 
is achieved within value ranges found in a 2D electron gas. Both device setups apply for arbitrary 
incoming electron polarization and operate at broad energy ranges within the incoming electron 
band. 

PACS numbers: 72.25.-b,85.75.-d,03.65.Vf 
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I. INTRODUCTION 



The Rashba and Dresselhaus spin-orbit (SO) interactions arise in materials which lack ei- 
ther structural or bulk inversion symmetry, respectively l|-|3|. These two kinds of interactions 
have recently been given a great deal of attention due to their potential role in thegeneration 
and manipulation of spin polarized currents, spin filters ^-[3|, spin accumulation |8|], and spin 
optics[9]. 

Recent proposals have been made for the construction of )erfect spin filters based on ac- 
tive Rashba spin orbit media[6] , ballistic spin interferometers 10| and the analysis of the per- 



sistent spin helixl? , . Further applications based on the interference concept include quan- 

n ^ n 

tum logic gates|l2|. bit controlled Stern-Gerlach devices[5] and tunable entanglement [13|. 
Here we readdress the problem of spin filtering by interferometry explicitly implementing the 
original concept in reference js, gI, proposing a quasi two dimensional device, to test the spin 
filtering concept through an electronic Mach Zehnder interferometer (MZl) within Rashba 

n n 

and Dresselhaus media|14j. |15| . Generally, spin injection is achieved by drawing currents 
from oriented ferromagnets, where the majority carriers have a preferred spin orientation. 
In order to change spin orientation it is then necessary to use strong magnetic fields and 
depend on spin relaxation times that limit manipulation speeds. The currents drawn by 
these methods still contain a large fraction of the undesired spin orientation depending on 
the magnetic field and the density of states. For this reason it is very desirable to achieve a 
larger fraction of oriented spin albeit not a pure spin current. Conceptually, the proposed 
device would be able to separate spin-orientation, as in a Stern-Gerlach device, but without 
magnetic field gradients, and both in the absence (asymmetric interferometer) and presence 
of a weak magnetic field (symmetric interferometer) coupled to the spin-orbit inter action [l^. 

The aim of this work is to generate spin polarized currents making use of a combination of 
weak magnetic fields (where the Zeeman term contribution is small), the intrinsic Dresselhaus 
spin-orbit interaction and the tunable structure inversion asymmetry (SIA) Rashba type in- 
teraction. The interference of electron precession and phase changes in the wave function 
render sufficient degrees of control to separate spin components into distinct output channels 
in a variety of experimental parameter ranges. Within this setup, we obtain precise condi- 
tions for spin filtering by numerically deriving the spin polarized electric currents generated 
by electrons drawn from an unpolarized reservoir as a function of temperature and chemical 
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potentials. The experimental conditions are derived using the Landauer-Buttiker transport 
formalism within a tight-binding model that incorporates realistic material parameters and 
the effects of voltage reservoirs. This description affords specific experimental parameters for 
building the filtering device regarding spatial dimensions, strength of material parameters 
and of the externally applied fields and biases, in order to achieve optimal spin polarized 
yields. 

It is iniportant to point out that the translation operator approach to this problem in 
reference jg, [l^ only provides a conceptual proof that in principle interferometry can yield 
spin selection. There, exact solutions can be found ignoring many practical drawbacks that 
should be evaluated for a successful proposal of a device, namely: The treatment presented 
here accounts for a finite input band of energies due to leads connected to a reservoir and 
to an output sink. It also accounts for the effect of temperature through the Fermi function 
at the reservoirs as long as the coherence length does not become shorter than the device 
size. Finally, finite voltage differences are contemplated so that an actual current can be 
driven through the device. The latter source of control of the input energies lends itself as 
an additional tuning parameter to optimize spin filtering, a tuning not accessible in previous 
more idealized proposals. 

We present two general device configurations: i) a symmetric arm configuration, that 
requires a weak magnetic field and renders an almost fiat energy dependence of the polar- 
ization conditions and ii) an asymmetric interferometer, that requires no external magnetic 
field and only Rashba couplings, with a more complex energy dependence that demands 
specific voltage ranges to filter. This proposal shows the practical possibilities, within a 2D 
electron gas, of the interferometric concept advanced theoretically by Hatano, Shirasaki and 



Nakamuraj6| 16| 



II. METHODS 

We consider a quasi two dimensional electron gas consisting of non interacting electrons 
subject to both Rashba and Dresselhaus spin orbit interactions. In addition, an external 
magnetic flux $Bj described by a vector potential A, threads the device transversally. Recent 
works have shown how to measure and control the Rashba parameter using gate voltages in 
two dimensional GaAs/AlGaAs electron gas 17H20| and also in other heteroj unctions, such 
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as InAs/AlSb and HgTe[21]. In general the Rashba parameter a = b{E), where (E) is the 
expectation value of the electric field at the 2DEG, and b depends on the inverse of both the 
effective mass and the material gap 



22|. Measurements of SO parameters have been made 
by either Shubnikov-de Haas or weak localization/antilocalization effects. While intrinsic 
Dresselhaus parameters cannot be changed by a gate potential, they can be tuned by strain 
effects|23|. 



One can address the two dimensional GaAs/AlGaAs electron gas by a single particle 
Hamiltonian including the previously described couplings by 

H = ^ + U - ^{U^ay - Uyan - f (n,a^ - U^a^ + (1) 

where U = {p + eA) 1 is the kinetic momentum, 1 being the 2x2 identity matrix (in 
spin space). Ordinary vectors are represented with an over-arrow while 2x2 matrices are 
written in bold face. The electron's charge and effective mass are denoted as — e and m*, 
U = f/ 1 is a substrate lattice periodic potential, a is a vector of Pauli matrices, and a and /3 
are material-dependent parameters characterizing the Rashba and Dresselhaus interactions, 
respectively {2^. The simple form for the Dresselhaus term is due to the assumption of strong 



confinement of the electron gas, so that kp << 7r/o?[24| where kp is the Fermi wavector of 
in plane electrons and d in the confinement length scale. The condition essentially states 
that the direction perpendicular to the gas can be averaged, reducing cubic contributions to 
the dominant linear terms described here. The Dresselhaus term depends on the orientation 
of the crystal axes, our geometry insures that all interferometer arms are described by the 
same Hamiltonian. The last term is the Zeeman energy, which we will ignore in the limit 
of small magnetic fields (a few fiux quanta through a 200 x 200/im^ area) where it is much 



smaller than the spin orbit energy 
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For aGaAs heterostructure a ~ 3.9xlO~^W m[26'|, /3 ~ 2.4xlO"^2eV m and h^/m*L 
1.7 X 10~^^eV m, assuming electron ballistic propagation of length ~ 1/im and an effective 
mass of m* = 0.067mo, with mo the free electron mass. The proposed device configuration is 
depicted in FigUl a Mach Zehnder Interferometer (MZI) implemented by voltage gating two 
interfering electron paths. GaAs heterostructures are known to have charge phase coherence 

n 

lengths of 20yum at 15mK while spin coherence lengths are larger: ~ lOO/im at 1.6K[5|. 
Thus, at 15mK both coherence lengths will be larger than a device of 20/xm size. 

We are interested in determining the polarized currents Id- at outgoing channels Di, with 
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FIG. 1: Electronic Mach Zenhder interferometer setup on a SO active GaAs/AlGaAs 2D electron 
gas and the corresponding tight-binding model superposed. The device is biased on the left at 
voltage +y/2 and the current is collected on the right at the two reservoirs biased at —Vj^. The 
electron beam splitters are implemented through two Schottky gates with reflection and transmis- 
sion coefficients p and r separating electrons between upper u and lower I tight binding arms. There 
is a magnetic flux <I> b perpendicular to the plane. The outgoing currents 1\ and I2 are collected at 
D\ and D^- 

i = 1,2 and to find the conditions for spin filtering at either output channel. We define 
the spin filter as one acting on any entering polarization and returning a polarized state along 
a definite axis. This approach will serve to build two polarized spin currents of opposite 
polarization. The relevant processes within the interferometer are described as follows (see 
figure [I]): Single electrons are assumed to be extracted from the left electron reservoir at 



voltage +V/2 
Schottky gate 



?he electrons then pass through the first beam splitter, implemented by a 
27| labeled (SGi)[28.] resulting in two intermediate output beams. These two 
beams interfere at the second Schottky gate (SG2), from which two final beams emerge at 
the output channels. 



As we consider both Rashba and Dresselhaus interactions, we need to derive general re- 



flection conditions at the Schottky gates. In reference |28|, this was done for the Rashba 
interaction, assuming that small enough spin-orbit strength would yield only a small diver- 
gence of the reflected spin states. In this paper we will take the case of 7r/4 reflections, that 



can be controlled experimentally, and lead to simple, spin-orbit independent reflection and 



transmission matrix elements 



14l |. If it were not the case and large deviations from such 



conditions arose the beam splitter would behave differently as the spin-orbit parameters 
changed and would complicate predicting an appropriate operating regime. 

A tight binding model will allow us to quantitatively parametrize for the two dimensional 
gas and spin-orbit intensities appropriate to a GaAs/AlGaAs junction. The tight-binding 
Hamiltonian for the interferometer SO active region is 



H 



TB 



E(4t 4 




(2) 



rX 



fX,r'X' 



where are the site energies, cl^ creates electrons at site r with spin A =t, i and t^'i is 
the transfer integral between sites/spins r*, A and r'. A'. The sum over the sites r draws 
the legs of the interferometer, like in Fig. [1] The states are denoted ) with 

^rx = 0, 1 depending on whether or not an electron with spin A occupies site f. The link 
between equations and ([2]) is established through the definition of the spinor components 
i/j^[r) = {fXlifj). 

The transfer integral between sites r and r' including both field and SO effects measures 
the phase aquired by the electrons when they are transported along the interferometer legs. 
It can be written in the "position-spin" basis {|rA)} as 



t^^/ = {fX\t exp 



|r'A') 



(3) 



where t is the coupling between sites and it is set to 0.156 eV (see Fig{T]), in accordance to 
the Fermi velocity in GaAs, (pQ = h/e is the flux quantum and 

1 



A{x,y) 



m 



+ 



1 777* 



(4) 
(5) 



accounts for the spin spatial-precession and phase changes due to the SO coupling and the 
magnetic field, that operates in between tight binding sites e.g. at the middle of the link 
which connects sites r and r'. We have chosen the simple staggered gauge for the magnetic 
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field vector potential. Equations (jSj) and 
described in terms of a minimal coupling 



5h are easily derived when the SO interaction is 



For the current incoming from the left electron reservoir and outgoing to the detectors, we 
have transmission and reflection coefficients r, p and r',p' respectively (see FigH]), that can 
be gate controlled. The gate is spin inactive as it constitutes a scalar potential barrier, for 
both the input and output channels, so that in the notation of Eq|3l t^^^J.^^ = r, t^^^^ o ~ P'> 
t^'^ n = P' ■, t^'^ n = t', t^'^ r, = t', t^'^ n = p' where r„o and fm are the first tight-binding 
sites of the arms and Tun and tin are the A^th sites of the corresponding arms. We choose 
throughout p = p' = i/ \/2 and r = r' = 1 / \/2. ^ ^ 

The coupling to the metallic reservoirs is described by a self-energy [35]. The self energy is 
derived by solving the Dyson equation for a chain of sites indexed by n where S„ = t ^^ t 
and S„ = for site energies chosen to be zero, for an infinite chain. The resulting 

quadratic equation can be solved and has real and imaginary parts so that 



- 2 '^IcadA^leadA' K'^) 

where we have introduced creation and annihilation electron states at each lead of the system 
i.e. in, Di, and D2 for the incoming lead and the two drains, where a metallic lead couples 
the system to the reservoirs. Three such couplings will be considered below, at the incoming 
beam and at the two output leads connecting to detectors. The transmission between a spin 
state A at the input and spin state r at the output Di, at energy E is given by 3^ 

T^^ = Tr [T^^G{E) G{Ey] , (7) 



D 



where the broadening, due to the coupling to the electron reservoirs is given by Tr 



— (Sr^)^ and the Green's function is computed as 

^(^) = E - - (S" + S^i + E^2) ' 
where the self energies are diagonal matrices that include all spin orientations at the corre- 
sponding leads jsG]. The current at each output lead Dj with spin component r is computed 
using 

= T E / ^d]{E) {r{E) - f^^{E)) dE. (9) 



The Fermi occupation at temperature T is given by f^{E) = l/(exp(-^j^^) + 1), where 
E'p = ±V/2. 



The numerical calculation would proceed in the following way: The incoming and outgo- 
ing one dimensional leads, are made up from site energies set to zero as a reference energy 
and nearest neighbor couplings t set to 0.156 eV. The self energy of such leads can be com- 
puted exactly 35|] leading to an energy contribution at the input and output sites of the 
interferometer. The Hamiltonian of the interferometer is then computed according to EqJH 
and the hopping matrix elements in Eq|3l The spin orbit parameters a and (3 are only 
different from zero within the interferometer. One can then combine the computed matrices 
for the Hamiltonian and the self energies, into the Green's function defined by Eq|8] which 
involves a matrix inversion and contains all the information on the external contacts of the 
interferometer. Such an expression is also a matrix and is used to compute the transmission 
between specific spin components through the interferometer by way of EqlTl The trans- 
mission is a number once the spin components have been fixed, and it must be computed 
for each of the energies within the incoming band. From this function, the spin polarized 
currents can be obtained. 

As a measure of the performance of the device we use the transmission asymmetry at a 
particular output channel, which is defined as 

An^ (0, a, /3) = Tl] + T^] - (T^^ + T^]), (10) 

considering that the states at the input lead are thermal mixed states at temperature T. In 
this case, the amplitudes do not interfere and the calculation requires summing probabilities. 
This measure is used due to the fact that, in contrast to the translation operator approach of 
previous proposals, we cannot obtain exact conditions for spin filtering when reservoirs are 
considered, although the arms of the device are one dimensional. Therefore, the conditions 
for operation must be sought numerically, with an adequate measure for performance. 



III. RESULTS 



Figure E] shows the spin filtering asymmetry at output lead Di as a function of the Rashba 
and Dresselhaus values for the material. The magnetic field has been chosen so that the 
spin asymmetry is maximized. The potential points for operation are identified in the figure 
within less than one flux quantum through the device. Such operation point comprises an 
asymmetry of more than 99% so the filter is highly efficient. Furthermore, the sensibility of 
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FIG. 2: Transmission asymmetry as a function of the Rashba and Dresselhaus spin-orbit interaction. 
The clearer shades represent positive asymmetry indicating up spin preference according to Eg 110} 
the values on the lines represent the magnitude of the asymmetry. Left panel: The symmetric 



interferometer (see inset), where t 
in GaAs/AlGaAs heterostructures 



le point of operation is indicated by a at values encountered 
2^ while X is a representative point for small asymmetry. The 
magnetic flux at the operation point is $ = 0.248(/i/e) a field of 10 Gauss through an interferometer 
of Ifirn^. Right panel: The asymmetric interferometer (see inset) where the point of operation is 
indicated by a (^. There is no applied magnetic flux at the operation point and the Dresselhaus 
parameter can be close to zero for an interferometer of 1/im^. 



the filtering asymmetry to changes in the parameters is weak, i.e. a 15% change in Rashba 
couphng only changes the asymmetry in 10% around the operating point. Thus, tuning 
around the operation point is by no means beyond experimental accuracy. 

Figure [2] shows the transmission at selected points in the a, (3 parameter space (see Fig|2]) 
as a function of the magnetic flux through the device from the input to the output channel 
Di. The top panel depicts an a,P choice with negligible filtering properties, where all 
transmission probabilities are appreciable so that only a weak asymmetry results. The 
bottom panel shows the transmission probabilities at the filtering operating point. The 
figure shows that one can choose the magnetic flux so as to obtain optimal filtering for 
either spin component (there are two optimal operating points in the range chosen for the 
SO parameters). The other components for the specific operation point are negligible for all 
fields shown. The signal at output channel D2 will show the opposite behavior as depicted 
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FIG. 3: a) Transmission probabilities as a function of magnetic flux for the point X in Figl2]left 
panel, off the filtering operating point, b) Transmission probabilities tJJ^ and T^_^ at lead Di at 
the operation point (see Fig l2] left panel), as a function of the magnetic flux. The transmission 
coefficients Tjj'_^ and T^^ are negligible at the operation point. The behavior of the lead D2 is 
exactly opposite to that of Di allowing for a separation of opposite components in a single device. 

by the dashed blue lines in the bottom panel [b) of Figj3] (exchanges fo^' ii); so we 
can simultaneously and separately draw both spin orientations. We note that there are no 
optimal operating points for either a or /3 = 0, for which only a small asymmetry results. 

An important feature of the proposed device is that, in the symmetric arm configuration, 
the filtering properties are robustly energy independent in the tight-binding model i.e. the 
properties remain in the full band width of the incoming channel[6]. Small interferometer 
arm asymmetries are nevertheless inevitable in a real setup, so as we will see below, an 
energy dependence is expected. Due to the insensitivity of the filtering to the input energy, 
the polarized current versus applied voltage is linear (see Eq. [9]) whenever ksT <^ eV i.e. 
the thermal energy is small compared to the potential difference applied to the device. 
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FIG. 4: Transmission probabilities T^^ and Tf^^ at lead D2 as a function of energy at the operation 
point in Figl2] right panel. The transmission coefficients T^^ and T^^ are negligible at the 
operation point. The input gate potential (Fermi energy) is centered at the peak of transmission 
for the up spin in D2, then the potential difference between the input and lead is increased as 
shown. 



The second configuration we analyze consists of an interferometer with a large asymmetry 
in arm lengths. In this configuration one can filter in the absence of a magnetic field using 
the difference in translational phase. As we have fixed the length difference, one does not 
achieve as high a filtering capability as in the symmetric arm configuration. This setup 
also compromises the energy dependence i.e. it is only for certain ranges of energy, within 
the input band, that the interferometer will work. Nevertheless, arm lengths could also be 
tuned to yield a higher asymmetry. Figure [2] right panel, shows the transmission probability 
asymmetry for the second interferometer. Again, one can tune the SO parameters in the 
same range as in the symmetric arm setup. Nevertheless, for the particular configuration 
chosen (a 3yum and a 1/im arm) the optimal point of operation occurs at zero Dresselhaus 



23|. 



coupling. The latter coupling could in principle be achieved by means of strain 

In Figure m we show the transmission probability as a function of the energy for the spin 
flipping processes (f— '"i and i— T-f), the other components being negligible. The transmission 
is no longer optimal in the whole energy range but peaks in certain energy intervals. If one 
tunes the input Fermi energy and then samples potential difference in the vicinity of the 
optimal transmission, for a particular output lead, one can obtain a substantial polarized 
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FIG. 5: Polarized current components at lead D2 as a function of the potential difference as indicated 
in FigjH The applied magnetic flux is zero. Above 0.1 V there is an increasing contribution from the 
opposite spin orientation that degrades polarization. The inset shows the proposed configuration 
that differs from that of Fig|T]by the position of the output leads. 

current within a certain bias range. Figure [5] shows the polarized current components for 
the situation depicted in FigHl The figure shows that one can achieve at least four times 
as much spin up polarized current as spin down current contribution within a 0.2 volt bias 
around an optimal transmission probability. 

IV. SUMMARY 

We have proposed a practical spin filtering device based on a Mach-Zehnder type spin in- 
terferometer built into a 2D GaAs/AlGaAs electron gas. Spin filtering operation is achieved 
by separating spin up and spin down components of an incoming thermal mixed state ex- 
tracted from a biased electrode. Spin filtering is possible by tuning Rashba and/or Dressel- 
haus spin-orbit contributions by way of an applied gate voltage or strain. The Landauer- 
Buttiker formulation, properly accounting for coupling to reservoirs at finite temperatures 
and voltage bias, is used to compute transmissions and polarized currents and demonstrate 
the robustness of the interference concept for current polarization. Two device modes are 
discussed that operate at complementary experimental conditions: i) A symmetric device 
that requires a weak magnetic field and tuning both Dresselhaus and Rashba interactions. 
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Such configuration displays a flat incoming electron energy dependent filtering, and ii) an 
asymmetric device that only requires tuning of Rashba couplings and operates at specific 
voltage biases. Such energy dependence makes the latter device a simpler setup but more 
prone to voltage and temperature sensitivities. Both configurations display the drawing of 
opposite spin polarized currents at the output leads at the optimal operation point. 

While the gate voltage control is well a established mechanism, strain has recently been 



tested for spin control by building mismatched layers in AlGaAs|37| generating spin splittings 
of the BIA type of around 0.007 meV, and 0.025 meV for the bulk Dresselhaus interaction. 
Strain can account for a 20% change in the BIA coupling in mismatched AlGaAsjsSl. This is 
within the requirements of the proposed device. An additional tunable parameter is a weak 
perpendicular external magnetic field that induces effects similar to tuning the interferometer 
arm lengths. The interferometric setup, proposed here, avoids using ferromagnetic material 
or strong magnetic fields (Zeeman term or spin torque mediated) to polarize electron spins. 
The device also operates at broad energy ranges for the incoming electron band, a fact which 
averts the need for difficult fine tuning of the required parameters. 
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